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Abstract 

Edge-based elements have proved useful in solving electromag- 
netic problems since they are divergenceless. Previous authors have 
presented several two-dimensional and three-dimensional elements. 
Herein, we present four types of elements which are suitable for mod- 
eling several types of three-dimensional geometries. Distorted brick 
and triangular prism elements are given in cartesian coordinates as 
well as the specialized cylindrical shell and pie-shaped prism elements 
which are suitable for problems best described in polar cylindrical 
coordinates. 
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where E is the electric field and G is a Green’s function. For the remainder 
of this report, we suppress the dependence of F on the unprimed coordinates 
since there will be no operations carried out in these coordinates. Further- 
more, we express (1) in terms of its cartesian components 

I = il x + yl y + zl z (3) 


where 

= C ' Jv E{x',y',z')dx' dy' dz (4) 

and ( € {x,y, z}. 

To evaluate (4), we apply the change of variables theorem to map the 
original brick domain into the unit cube as shown in figure lb. 


k~C- f / / F(x'(u,v,w),y'(u , 
*/ o »/o J 0 

which is accomplished by allowing 

x’(u,v,w) = 

y'(u, v, w) = 
z'(u,v,w) = 


u>)» •*'(“, u »w))|.7(u, 

v,w)\dudvdw (5) 

8 

J2 N i(u, V, w) x\ 

(6) 

1=1 

Y 2Ni(u,v,w ) y'i 

(7) 

t = 1 

Ni ( u > v > w ) z i 

i=i 

(8) 


with Zj) denoting the ith node coordinates as shown in figure la. The 

shape functions are given by 


JV i(u,v,w) 
N 2 (u,v,w) 
AV'u, v, w) 
N 4 (u, v, w ) 


(1 — u)(l — u)(l — w) 
u(l — v)(l — w) 
uv(l — w) 

(1 — u)v(l — u>) 


Ns(u, V, to) 
N 6 (u,v,w) 
N 7 (u, V, to) 
N%(u,v,w) 


and the Jacobian is of the form 


(1 — u)(l — v)w 
u(l — v)w 
uvw 

(1 — u)vw 

(9) 


dx' dy' dz' dx' dy' dz' dy' dx' dz' 

du dv dw du dw dv du dv dw 

d£dx^_d£ dE_d£dy^_ _dz^dx^dy^_ 

+ du dw dv du dv dw du dw dv 


( 10 ) 
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1 Introduction 

Edge-based expansion functions have become popular in finite element and 
moment method techniques applied to the solution of electromagnetic scat- 
tering problems. This is due to the fact that the expansion functions are 
divergenceless within the element domain while any fictitious charges formed 
at the edges of the element will be canceled by another edge resulting in a 
volume which is divergenceless. In contrast, nodal expansion functions are 
not divergenceless and are thus less accurate. 

Tanner and Peterson [1] presented triangular and tetrahedral divergence- 
less vector elements while Jin and Volakis [2] applied the corresponding rect- 
angular brick elements to slot scattering problems. Kameari [3] presented 
distorted linear and quadratic brick elements in the context of transient 
quasi-magnetostatic problems. 

In this report we derive edge-based expansion functions for a variety of 
three-dimensional elements. These expansion functions are characterized bv 
unit value along one edge of the element while vanishing along all the re- 
maining edges which are parallel to the unit edge. Distorted bricks, distorted 
triangular prisms, cylindrical shells and pie-shaped prisms will be presented. 
The first two of these are suited for completely general inhomogeneous do- 
mains whereas the last two are more appropriate for domains terminated at 
cylindrical boundaries. 

2 Distorted Brick Element 

The first element considered is the distorted brick which is shown in figure la. 
The strategy used herein is to first employ the change of variable theorem to 
convert the integration domain from the original distorted brick into a cube. 
Once this is accomplished, the vector-valued function may be expanded in 
terms of edge-based expansion functions as was done by Jin and Volakis[2]. 
Typically, one wishes to integrate some vector- valued function 

I{x,y,z) = J v F(x',y',z',x,y,z)dx' dy' dz' (1) 

over the domain of the element where V is the volume of the brick and the 
vector valued integrand is typically of the form 

F{x\ y\ z) = E{x', y', z) G{x', y', z\ x, y, z) (2) 
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where E is the electric field and G is a Green’s function. For the remainder 
of this report, we suppress the dependence of F on the unprimed coordinates 
since there will be no operations carried out in these coordinates. Further- 
more, we express (1) in terms of its cartesian components 

/ = xl x + yl y + zl t (3) 


where 

h = C • J v E(x',y',z')dx' dy' dz' (4) 

and C € {x,y,z}. 

To evaluate (4), we apply the change of variables theorem to map the 
original brick domain into the unit cube as shown in figure lb. 

/c = C / / / F(x'(u,v,w),y'{u,v,w),z'(u,v,w))\J(u,v y w)\dudvdw (5) 

0 wO J 0 

which is accomplished by allowing 
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x'(u , V, w) 

= ^iV,(u,y,u;)z- 
1=1 

(6) 

y'(u,v,w ) 

= 53 N <( u ' v ' y'i 

(7) 


t=i 
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z'(u, V, w) 

= J2 N <( u ’ v ’ w ) z 'i 

»=l 

(8) 


with (x,, y iy z t ) denoting the ith node coordinates as shown in figure la. The 
shape functions are given by 


Nx(u, v, w) 
jV 2 (u, v, w) 
N 3 (u,v,w) 
iV 4 (u, v y u;) 


(1 - u)(l - u)(l - w) 

u(l — u)(l — w) 
uv(l — w) 

(1 — u)u(l — w) 


N s (u,v y w) 
N 6 ( u,v,w) 
Nt(u,v,w) 
N&(u, v,tv) 


and the Jacobian is of the form 


(1 — u)(l — v)w 
u( 1 — v)w 
uvw 

(1 — u)vw 

(9) 


dx’ dy' dz' dx' dy' dz' dy' dx' dz' 
du dv dw du dw dv du dv dw 

dzl_d£dy^_ _ d£dx^_dy^ 
du dw dv du dv dw du dw dv 
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( 11 ) 


The partial derivatives are computed from 


d£ = A 
^ & 


dN if .. 
— {u,v,w) C, 


with C € { x,y,z } and v 6 { u,v,w }. From (9) the partial derivatives of the 
shape functions are found to be 


dNi(u, v s w) 

du 

= — (1 — w) (1 — tv) 

dN S {UsV,w) 

du 

= 

— (1 — v) w 

5iV 2 (u. v, w) 

= (l-v)(l-tl>) 

dN 6 (u,v, tv) 

du 



(1 — v) w 

du 



dNz{u,v,w) 

du 

= v (1 — w) 

dN 7 {u , v, w) 

du 

= 

vw 

dN 4 (u , v , u>) 

du 

= -v(l — tv) 

dNs(u , v , w) 

du 

= 

—vw 

dN x {u , v, tv) 

dv 
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dN 5 {u,v,w) 

dv 

= 

— (1 — u) w 

dN 2 (u, v, tv) 

dv 

= -u (1 - w) 

dN 6 {u,v,w) 

dv 

= 

— UU) 

dN 3 (u, v. tv ) 

= u (1 — w) 

dN 7 (u,v,w) 

dv 

= 

UIV 

dN 4 (u, v, w) 

dv 

= (1 — u) (1 — w) 

dN 6 (u,v,w) 

dv 

= 

(1 — u ) w 

dN\{u,v,w) 

dw 

= -(l-u)(l-v) 

dNs(u,v,w) 

dw 


(l-u)(l-f) 

<9JV 2 (u, u, tv) 

dw 

= — u (1 — u) 

dN e (u,v,w) 

dw 

= 

u (1 — v) 

dN 3 (u,v,w) 

duo 

= —UV 

dN 7 (u,v,w) 

dw 

= 

UV 

dN 4 (u,v,w) 

dw 

= - (1 - u)v 

dN s {u,v,w) 

dw 


(1 — u ) V 
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Now that we have transformed the integration domain from the original 
distorted brick into the unit cube, we may expand the vector-valued integral 
in terms of an edge-based expansion of the form 


12 

F{x'(u,v,w),y\u,v,w),z'{u,v,w)) = N*(u, u, w)ipj (12) 

j=i 


where 


n;(u, 

V, 

w ) = 

(1- 

v)(\- 

w)t u 

iV 7 '(u, 


w) = 

u(l - 

™)t 58 

NJ(u, 

V, 

, w) = 

v(l - 

- w)U 3 


N!(u, 

v, 

w) = 

uwt 67 


Ni(u, 

V, 

, w) — 

(1- 

v)wt 56 


Ni |(u, 

U, 

w) = 

(1 — U 

0(1 - v )i 

Nl(u, 

V, 

w) — 

vwt 87 


N io( u , 


w) = 

u(l — 

V)t 26 

#}(«. 

v, 

w) = 

(1- 


w)t u 


v, 

w) = 

(1 — U 

00 
< -vj 

Ni(u, 

V, 

w) = 

(1- 

u)wt 23 


N^(u, 

U, 

w) = 

uvt 3 7 
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in which xpj are the unknown expansion coefficients and t XJ denote the ele- 
ment’s unit edge vectors given by 



with fj being the position vectors in the original domain. 
We may now express the cartesian components of (3) as 


iV*(u,u,u;)|J(u,u,iij)| du dv dw 


(14) 


(15) 


recognizing that ipj represents the average component of the vector-valued 
function associated with the jth edge. 

It should be noted that the brick cannot be too distorted. This is due to 
the change of variables theorem. In order to assure a one-to-one mapping for 
(x, y, z) to (u, v, w), the sign of the Jacobian may not change. This condition 
will preclude elements that axe twisted or extremely narrow and diamond 
shaped. The typical element should therefore not deviate much from the 
rectangular brick. 
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3 Distorted Triangular Prism 

A triangular prism element such as the one shown in figure 2a may be de- 
veloped in exactly the same manner as the brick. In this case, we employ 
a mapping from the original volume V to the volume shown in figure 2b 
followed by an edge-based expansion of the vector-valued integrand. The 
cartesian component integral may be expressed as 


/, = (• /V 1 / 1 * F{x'{u,v,w),y'{u,v,w),z'(u,v,w))\J{u, 

Jo Jo Jo 

u, tu)| dvdu dw 
(16) 

where the coordinate mapping is given by 


6 

x'(u,v,w) = ^TiV,(u,t>,u;)x' 
1 = 1 

(17) 

6 

y'(u,v,w) = ^2Ni(u,v,w)y' t 

1 = 1 

(18) 

6 

z'(u,v,w ) = 53 iV ‘( u ’ w ) z i 

1=1 

(19) 

with the node coordinates (x,, y t , z { ) shown in figure 2a and the shape func- 


tions are given by 


Ni(u,v,w) = (1 - tt - u)(l - w) N 4 (u,v,w ) = (1 - u - v)w 

N 2 (u,v,w ) = it(l-u>) N 5 (u,v,w) = uw 

N 3 (u,v,w ) = v(l-u>) N 6 (u,v,w) = vw 


The Jacobian is again given by (10) with 


( 20 ) 


and the non-zero shape partial derivatives of iV, are 


dNi 

= -(l 


dN^ 

dNi 

= -(l- 


du 

— w) - 

dv 

dw 

u — v] 

dN 2 

= (1- 


dN 3 

dN 2 


dNs 

du 

w) = 

dv 

dw 

= -u = 

dw 

dN 4 


_ dN 4 


dN 3 


dN 6 

du 

= — w 

dv 


dw 


dw 

dNs 

= w = 

dNe 





du 

dv 






dN* 

dw 


( 00 ) 


The vector valued function in (16) is expanded in terms of edge- based func- 
tions 

9 

F{x'(u, v , w),y\u, v, w), z'(u, v, u;)) = N‘(u, v, w) xj>j (23) 

i=i 

where the vector expansion functions are 


jV* (ti, v , w) 
N% (u, v, w) 
N£(u,v,w) 
N%(u, v, w) 
N%(u,v,w) 
N‘(u,v,w) 
N‘(u,v, w) 
JV|(u,u,u>) 
JV‘(u,v,w) 


[(6 + £2)^12 + ^2*13](1 — U >) 

— \/2[^3<12 — ^ 2 ^ 13 ] (1 — w ) 

— [^3^12 + (£l + ^ 3 )^ 13 ] ( 1 ~ w) 
[(£1 + ^2)^45 + f 2 * 4 e]w 

— >/2[^3^45 + ^2^46]^ 

~ [£3*45 + (4l + 6 j)£m]w 
u(l - u)t 28 
(1 - u)(l - v)ti 4 
(1 - U)vt 36 


(24) 


where the area coordinates in (24) are given for the new integration domain 

£n = 2^Sn(s n <2 n )(a n b n )(s n Cjj) 

■®n = 2 (®*» 4" b n -f- Cn) (25) 

«€ {1,2,3} 
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and 


a x = 

\/(« - l) 2 + V 2 

a 2 = 

\Ju 2 + (v — l) 2 

a 3 = \/u 2 + v 2 

bx = 

Ju 2 + (v — l) 2 

i>2 = 

yju 2 + V 2 

b 3 = y / (u-l) 2 + f 2 

C\ = 

V2 

c 2 = 

1 

c 3 = 1 





(26) 

We are 

now able to express 

(16) as 




k = X>; J J J N‘(u,v,w)\J(u,v,w)\dvdudw (27) 


realizing that the same limited distortion constraint exists for the triangular 
prism as was the case for the brick. 


4 Cylindrical Shell Edge-Based Elements 

Although the brick and triangular prism elements given previously are quite 
general, we find it advantageous to develop an element which is designed 
for cylindrical problems. The resulting formulation will be quite compact 
compared to the brick for example since a specific geometry is assumed. This 
element is the polar cylindrical coordinate equivalent of the brick introduced 
by Jin and Volakis [2] and is shown in figure 3. 

For cylindricad problems, we wish to evaluate the integral 

/= rt*r F{ P \ t', z') P ' d P ' di,'dz' ( 28 ) 

J z\ J<t > i Jp i 

where the vector- valued function F(p', is analogous to (2). We do not 

employ any change of variable for this element so we may directly express 
the vector-valued function in terms of an edge-based expansion 

F{p\d,',z') = Y,N](p',d>',z')d,, (29) 

J = 1 


The shape functions axe given by 




1 

Ci. 


V - Zf> 

pI-p\\ 


Z\ ~ 25. 
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A 'W 

Nt{p\4>\z') 

Nl{ P \4>\z') 

Nt{p\<t>',z') 
Nl{p',4>\z') 
N e 7 {p\4>',z') 
Ni{p\4>\z) 
Rl{p',<tf>z') 
N‘ 0 (p',4>\z) 
Nn(p\ <t>',z') 
N' n (p',<t>\z') 


1 - 


|02 — <t>'\ 


102 — 0i I 



1 

1 

b* 

<jy 


m z l “ z 5 J 


V 2 - pi 


pi- Pi] 


z - z 5 


1 - 


Z\ - Zs 

10' - 01 1 
102 ~ 01 1 


0 



ID 

1 

"V» 

1 


m z l “ z 5 J 


P\ ~ P 


.P\~ P\ 



l 

u 

1 


.*1 - z s\ 


0 


1 ~ 


102 - 0'| 


1 02 ~ 01 1 J 



r— 

1 

"n 

1 


l^s-zij 


y 2 - p\ 

.p\- p\. 


Z_-_ Zj_ 
Zs - Zx 


0 


| 0 ' - 01 1 2 - Zi 

|02 — 01 1. _Z 5 — Z\ _ 


\P4-P'} 

[|02-0 , |1 

1 

.102 ~ 01 1. 

1 

Cl 

I 

■|0'-0lf 

H 

1 

■«r 

.102 — 01 1. 

1 

■»» 

1 

»— * 

10' — 01 1 

LP< ~ Pi J 

|02 — 01 1. 

P' ~ P 1 

" 1 02 0 7 1 " 

. P 4 ~ Pi. 

.1^2 “ 


0 


(30) 


where (p t , 0,, z,) denote the node coordinates as shown in figure 3 and (28) 
is given by 
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5 Pie-Shaped Prism 

The final element considered is the complement to the cylindrical shell ele- 
ment, the pie-shaped element which is shown in figure 4. The development 
of this element closely parallels the previous element. 

The vector-valued integral is given by 


/= rr r ns, *'>*') p'dpfd+'dz' 

J z\ J 4>2 J p i 

The vector-valued function is expanded in terms of edge vectors 


(32) 


j=i 


and the edge-based expansion functions are 

N{{p\<t>',z) = 

NZ(p',4>',z') = 

NZ(p\<t>',z') = 

NZ(p',4>',z') = 

Nl{p\4>\z') = 

Nl{p,<i>',z) = 

iV T V,*V) * = 


1 

w- 

4>i\ 

Za 

— 

2' 

I — 

. 

1 

CO 

<t>i\. 

Za 

— 

Z\ 

\p'- 

- Pi 

Z 4 - 

z'* 



.Pa ■ 

~Pl. 

Za - 

Z\ 

Y 


1 

\<f>3 ~ 

■4f 11 

za 

— 

z ' " 

1 “ 

\<h ~ 

<t>1\. 

Za 

— 

Z\ 


\£z 

<t>1\ 

z 

— 

Z\ 

I — 

\<f>3 ~ 

$1\. 

z A 

— 

Z\ 


-Pi 

’z 

Z\ 

6 

.P 1 

-Pi. 

za ~ 

■ Z X 

Y 

1 

|#3 ~ 

-f\] 

\i 

— 

Z\ 

1 — 

1^3 " 

- <fo|_ 

L*.- 

Z 1 

V 

-Pi 

|<Aa 


il 

y 

.p 1 

~Pi. 

.1^3 

- h\. 


\(f 

~ Pi 

1*' 

- <f>i\ 

A 

T 

.pi 

~Pi. 


- $i\. 

<4 


Nl{p\4>',z') = 


Pi- P 


Pi — Pi 


(33) 


(34) 
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where (/?,. c>,, -,) are the nodal coordinates shown in figure 4. With this 
expansion. (32) may be written 



6 Conclusions 

Four different linear volume elements were presented: distorted brick, dis- 
torted triangular prism, cylindrical shell, and pie-shaped prism. All of these 
elements are used to evaluate a vector- valued integral over the volume of the 
element by employing edge-based vector expansion functions. These func- 
tions were explicitly given for all elements. In addition, the shape functions 
and Jacobian associated with the mapping of the arbitrary brick or triangular 
prism to a known prism are presented. Since these elements are divergence- 
less, we expect them to useful for simulating various geometries via the finite 
element method or the Method of Moments. 
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